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i.  INTRODUCTION 


High  power  lasers  which  operate  in  the  infrared  (IR)  spectral 
region  are  used  in  advanced  systems  now  being  researched  and  developed 
by  the  Army.  There  is  a critical  need  for  efficient  components,  e.g., 
frequency  shifters,  spectral  filters,  modulators,  and  rotators  to 
realize  the  full  potential  of  these  lasers.  Conventional  modulator 
materials  such  as  KDP  operate  satisfactorily  in  the  visible  wavelengths 
but  are  ineffective  at  IR  wavelengths  for  use  in  high  power  IR  systems. 
The  term  modulator  is  used  here  in  the  broad  sense  to  mean  a device 
by  which  the  beam  can  be  manipulated  or  controlled.  Thus  a rotator  is 
a special  type  of  modulator.  Among  the  most  promising  conventional 
modulator  materials  in  the  IR  are  GaAs  and  CdTe  which  have  been  investi- 
gated at  10.6  u.  These  materials  have  shortcomings  for  amplitude 
modulation  and  large  index  phase  modulation  applications.  Due  to 
insufficient  linearity,  small  linear  and  angular  apertures,  and  high 
absorption,  they  are  unsuitable  for  use  in  powerful  CO^  laser  systems 

to  which  the  Army  is  committed.  These  applications  of  materials  utilize 
their  bulk  properties. 


II.  THEORETICAL  CONSI  DC  RATIONS 

A radically  different  approach  to  the  problems  of  designing 
low  loss  modulators  for  the  IR  is  shown  here.  This  approach  has  high 
merit.  In  this  technique  the  Faraday  rotation  of  Zeeman  split  energy 
levels  in  the  material,  due  either  to  doping  or  the  molecular  structure 
of  the  material,  is  used.  To  see  how  this  can  be  done,  the  complex 
index  of  refraction  n which  describes  optical  properties  of  a material 
is  considered.  The  real  part  n of  ft  is  associated  with  the  dispersion; 
the  imaginary  part  K is  a function  of  the  absorption  a 

K - £ a (1) 


where  cu  is  the  angular  frequency  of  light  and  c is  the  speed  of  light 
in  a vacuum.  Thus,  for  right  circularly  polarized  light  (+)  or  left 
circularly  polarized  light  (-) 


ft , * n,  - i K+ 


(2) 


Faraday  rotation  $,  the  angle  of  rotation  of  the  plane  of  polari- 
zation of  plane  polarized  light,  is  given  by 
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Figure  1 illustrates  how  magneto-optical  properties  of  bound  state 
transitions  result  in  a splitting  of  the  absorption  line,  which  gives 
rise  to  Faraday  rotation. 


Electric  dipole  transitions  which  change  the  magnetic  quantum  num- 
ber M by  +1  or  -1  correspond  to  right-  or  left-circulat  polarization, 

J 

respectively,  as  represented  in  Figure  1.  The  upper  state  degeneracy, 
removed  by  an  external  magnetic  field,  B,  is  indicated  together  with 
the  resulting  absorption-line  components,  which  are  shown  as  a function 
of  photon  energy,  E.  Since  the  absorption  a(E)  and  angle  of  rotation 
9(E)  of  the  plane  of  polarization  are  functions  of  the  imaginary  and 
real  part  of  the  complex  index  of  refraction  respectively,  6(E)  can  be 
calculated  from  oc(E)  by  use  of  the  dispersion  relation  [1] 


6+(E) 


g2  r a± (E* )dE' 

2v  JQ  E’(E2-E’2) 


(4) 


6(E)  = 0+(E)  - 0_(E)  has  been  calculated  from  assumed  shapes  of 
a+(E)  by  numerically  integrating  Equation  (4).  These  computed  values 

are  compared  with  the  measured  values  [2]  of  Faraday  rotation  due  to 
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Fg  -*■  transitions  in  CaF^iSm  . Excellent  agreement  between  cal- 
culated and  measured  rotation  spectra  was  obtained,  as  is  shown  in 
Figure  2 with  Gaussian-shape  absorption  components  a+(E)  and  oc_(E). 

Lorentzian  shapes  of  the  same  width  were  also  tried  and  did  not  fit  as 
well.  The  Faraday-rotation  spectrum  in  Figure  2(b)  was  calculated  by 
normalizing  a Gaussian  line  shape  to  the  measured  peak  absorption 
coefficient.  The  absorption  coefficient  y+  is  plotted  with  a dotted 

line  in  Figure  2(b)  and  is  related  to  a+  by  p+  = na+,  where  n is  the 

bulk  index  of  refraction,  taken  here  to  be  1.45.  Although  a closer 
fit  between  calculated  and  measured  curves  could  have  been  obtained 
by  adjusting  line-shape  parameters,  agreement  was  considered  already 
adequate  to  establish  the  validity  of  the  computational  procedure  used. 


The  rotation  9(E),  computed  for  two  different  separations  of 
a+(E)  and  a_(E),  is  shown  in  Figures  3a  and  3b  for  Gaussian-shape 

components  that  have  equal  amplitudes  and  half  widths.  At  a position 
Eq,  9(E)  decreases  less  rapidly  than  a(E)  as  the  separation  between 

the  two  components  increases.  The  computed  ratio  a(Eo)/9(EQ),  plotted 

as  a function  of  component  separation  6,  in  units  of  half  width,  is 
shown  in  Figure  3c  for  the  spectral  region  between  a+(E)  and  a (E) , 

where  o(E)  - a,(E)  + a (E).  It  can  be  seen  that  for  6 > 5,  a(E  )/9(E  ) 
+ — — o o 

is  very  small  compared  to  its  initial  value. 
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Figure  2.  Measured  rotation  spectrum  (a)  of  ^F.  - transition 
2+  U 1 

In  CaF^^m  at  4.3  k,  B ■ 60  kG.  Rotation  spectrum  (solid 

line)  calculated  (b)  from  Gausslan-shape  absorption  components 
(dotted  lines)  normalized  to  the  measured  peak  absorption 
coefficient. 
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Figure  3.  Rotation  spectra  represented  by  the  solid  line  in  (a) 
and  (b)  computed  for  different  separations  6 of  a+(E)  and 

a (£) , represented  by  dotted  lines  which  have  Gaussian  shapes. 

Numbers  shown  along  the  ordinate  must  be  multiplied  by  the 
maximum  values  of  a+  to  obtain  rotation.  (c)  Computed  ratio 

o(Eq)/0(Eq)  as  a function  of  component  separation. 


At  a separation  of  approximately  5 half  widths,  a wave  with 
frequency  w could  be  expected  to  propagate  through  the  material  with 

a very  low  absorption/rotation  ratio.  In  addition,  no  ellipticity  $ 
would  be  induced  by  the  material  into  the  wave  that  would  destroy  its 
linear  polarization.  This  result  follows  from  the  definition  [3] 
p = l/4[a+(E)  - a_  (E) ] , which  is  equal  to  zero  at  Eq  in  this  case, 

because  a. (E  ) = a (E  ) ~ 0.  The  absence  of  induced  ellipticity  is  a 
+ o - o 

desirable  feature  for  optical  isolators,  because  complete  isolation 
would  not  be  possible  if  the  wave  were  elliptically  polarized.  For 
an  effective  g of  9.6,  5 > 5 is  obtained  with  a few  kilogauss.  Effec- 
tive g values  as  large  as  9.6  occur  in  typical  bulk  transparent 

2+ 

materials  such  as  CaF^hy  in  the  visible  spectral  region  and  FeCo^ 

in  the  far  IR  [4,5];  thus  low-absorption  magneto-optical  devices 
appear  feasible  if  the  laser  frequency  is  matched  to 

The  Verdet  constant  A calculated  from  a reported  [5]  linewidth  of 

-1  2+ 

0.024  cm  and  an  effective  g of  9.6  in  CaF^iDy  , is  compared  in 

Table  1 to  conventional  rotators  which  utilize  a bulk  property  of  a 
material. 

A,  defined  as 


A 


e 

(XB) 


(5) 


where  X is  the  specimen  thickness,  was  evaluated  from  Figure  3(a),  which 

represents  a splitting  of  approximately  2.5  linewidths  or  0.06  cm  ^ in 
the  present  example. 

It  follows  from  Equation  (5)  and  Figure  3(a)  that 

(0.042)  ^ (60)  . (6) 

IT 

The  required  magnetic  field  B = 135  G was  found  from  the  line-splitting 

rate  of  0.45  cm  ^/kG.  It  is  assumed  that  (a)  = 10  cm  ^ and  from 

± max 

Equation  (3)  it  follows  that  A = 10  min/cm-G. 

The  calculated  insertion  loss  follows  from  the  Gaussian  line 
shape  of  the  components  and  the  maximum  value  of  a+(E).  Thus 


TABLE  1.  MAGNETO-OPTICAL  CONSTANTS  OF  MATERIALS 


T 

(°K) 

Verdet  Constant 
(min  G ^•cm 

A 

(microns) 

Insertion 

Loss 

(dB) 

TAG  [1] 

300 

1.27 

0.514 

0.4 

Corning  8363  [3] 

300 

0.1 

0.61 

1.0 

InSb 

77 

0.66 

10.6 

2.0 

2+ 

CaF^rDy  (Calculated) 

10 

10.0 

2.360 

0.0 

-5^  ' 10_^^  -1 

a (E  ) = (a  ) e [cm  ] at  6 = 5 . (7) 

± o ± max 


Only  absorption  due  to  transitions  between  bound  states  is  included 
in  this  calculation, 

A low  loss  IR  rotator  has  been  designed  based  on  the  preceding 
analysis  applied  to  data  of  Blum  et  al.  for  NO  [6], 


III.  DEVICE  APPLICATION 

Laser  radiation  of  wavelength  5.26  p is  directed  from  left 
to  right  in  Figure  4 through  a polarizer  onto  the  rotator.  An  absorp- 
tion line  in  the  rotator  is  split  into  two  components  by  the  applica- 
tion of  an  external  magnetic  field  from  the  superconducting  solenoid. 
The  strength  of  the  magnetic  field  is  adjusted  so  that  the  light  of 
frequency  passes  through  the  rotator  at  a spectral  position  between 

the  split  absorption  components.  Thus,  very  little  absorption  will 
take  place  and  the  plane  of  polarization  will  be  rotated  without 
introducing  ellipticity  into  the  beam  polarization.  The  length  of  the 
rotator  is  chosen  to  cause  a 45°  rotation.  An  adjustable  heater  and 
a thermocouple  are  provided  to  temperature  tune  and  stabilize  the 
rotator.  After  leaving  the  rotator,  the  light  beam  passes  through 
the  analyzer  which  is  orientated  with  its  plane  of  polarization  coin- 
cident with  that  of  the  light  beam.  Thus,  the  light  beam  will  pass 
through  left  to  right  but  not  in  reverse,  because  light  coming  from 
the  right  in  Figure  4 will  be  polarized  by  the  analyzer  at  an  angle 
of  45°.  After  passing  through  the  rotator,  an  additional  45° 
rotation  will  occur;  thus  a light  beam  entering  from  the  right  will 
be  polarized  90°  to  the  polarizer  and  therefore  will  not  be  trans- 
mitted through  it.  The  result  is  to  isolate  optically  everything  to 
the  left  of  polarizer. 


For  NO,  using  values  reported  by  Blum  et  al.  [6]  (left  side  of 
Table  2),  the  design  parameters  for  an  isolator  are  determined  (right 
side  of  Table  2). 


TABLE  2.  DESIGN  PARAMETERS  FOR  NO  ROTATOR 


Blum  et  al.  [6] 

Calculated  Design 
Parameters 

B 

30  kg 

100  kg 

Separation  of 
components  in 
halfwidths 

2.0 

5 

Rotation  per 
unit  length 

5 . 5 0 / cm 

2.4°/cm 

Cell  length 
for  45°  rotation 

18.8  cm 

Pressure 

3 torr 

3 torr 

Rotator 

absorption 

A1_-10  -1 
vlO  cm 

Correction  for  a finite  laser  linewidth  was  found  to  be  negligible 
for  a width  of  5,  the  maximum  value  assumed. 


IV.  DIAGNOSTIC  APPLICATIONS 


To  determine  their  suitability  as  modulators,  it  is  necessary 
to  investigate  basic  magneto-optical  properties  of  bound  state  transi- 
tions within  the  materials.  The  technique  [7,8]  developed  allows  the 
determination  of  such  parameters  as  g values  from  the  Faraday  rotation 
spectrum  which  previously  had  been  obtained  from  much  more  difficult 
high  resolution  absorption  measurements  [9]. 


Figure  5 shows  that  absorption  components  due  to  transitions 
7 5 5 

from  F to  D. , and  D , levels  were  not  resolved  with  the  minimum 

° 1 _1  -1  5 

usable  spectral  slit  width  of  0.6  cm  . Faraday  rotation  allows  D+^ 

level  splitting  to  be  determined  accurately  (Figure  6).  The  Faraday 
rotation  spectrum  was  resolved  into  two  dispersion  curves,  n+  and  n . 


1 
I j 
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gure  6.  Faraday  rotation  of  Fq  -*■  transition  in  CaF2:Sm 

£ is  the  separation  between  and  D_^  components  due  to 

an  external  60  kG  field. 


As  indicated  by  the  overall  shape  of  the  Faraday  rotation  spectrum, 
part  of  the  n+  curve  which  has  a positive  slope  in  the  region  G > 0 

is  unperturbed  by  the  n curve.  Because  it  is  known  from  the  disper- 
sion relation  that  the  n+  and  n curves  are  symmetrical,  the  shape  of 
the  n+  curve  which  has  positive  slope  in  the  region  0 < 0 is  also 
known.  The  remainder  of  the  n+  curve  was  determined  by  fitting  these 
two  known  parts  together  such  that  the  n+  curve  and  also  the  difference 
between  the  Faraday  rotation  spectrum  and  n+  were  symmetrical  about 
the  6-0  axis.  This  difference  is  the  n_  curve.  It  can  be  shown  from 

the  Kramers-Kronig  relation  that  the  peak  of  an  absorption  component 
occurs  at  the  spectral  position  where  the  corresponding  dispersion 
curve  crosses  the  axis  [10].  Subtracting  the  zero  point  of  one  dis- 
persion curve  from  the  other  in  Figure  6,  it  follows  that  the  splitting 

£ between  and  levels  is  2.2 ±0.3  cm  \ It  is  found  that 

g “ 0.39  ± 0.03,  from  the  Zeeman  shift  of  each  component,  given  by 


— S-r  gH  = 4.67  * 10"5  gH  - f , 

4-irmc 

where  £ is  in  cm  * and  H is  in  G.  This  result  is  very  close  to  the 

7 2+ 

theoretical  value  g = 0.40  for  the  state  of  a free  Sm  ion, 
calculated  from  the  Lande  formula 


„ J(J+1)  + S (S+l)  - LCL+1) 

8 2J (J+l) 

where  S - 3,  L = 7,  and  J » 4,  assuming  pure  LS  coupling. 

Close  agreement  between  the  g value  determined  from  Faraday 
rotation  and  that  calculated  for  a free  ion  supports  the  assumption  of 

2+ 

Kaiser  et  al.  [11]  that  these  transitions  in  CaF^cSm  occur  within 
highly  shielded  shells  (e.g.,  f shells). 

Results  from  Faraday  rotation  are  consistent  with  high-resolution 
absorption  data  reported  by  Zakharchenya  and  Ryskin  [9],  although  a 
somewhat  lower  g = 0.31  was  found  by  Margerie  from  magnetic  dichroism 
studies  [12].  This  technique  has  been  extended  to  treat  more  compli- 
cated energy  level  schemes  such  as  the  F-center  in  BaF^. 
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This  extension  of  the  computational  method  used  here  was  applied 
to  the  accepted  energy  level  scheme  for  an  F-center  in  a alkali  halide 
lattice  to  determine  the  type  of  defect  responsible  for  absorption 
bands  in  the  more  complex  alkaline  earth  fluorides,  e.g.,  BaF2»  CaF2» 

and  SrF2«  Experimentally  determined  Faraday  rotation  spectra  and  those 

computed  from  the  model  agree  closely  in  the  overall  features  for  low 
temperature  X-ray  Irradiated  BaF2«  As  before,  Gaussian  curves  gave  the 

best  fit.  With  the  aid  of  this  analysis,  it  was  found  that- the  accepted 
energy  level  scheme  for  an  F-center  in  a lattice  with  0^  symmetry  is 

adequate  to  explain  the  more  complex  case  of  a single  electron  trapped 

in  a lattice  with  T,  symmetry. 

a 

V.  CONCLUSIONS 

Calculated  and  experimentally  determined  plane-of-polarization 
curves  presented  demonstrate  that  magneto-optical  properties  of  transi- 
tions between  degenerate  energy  levels  may  be  used  to  design  low  loss 
components  for  laser  systems  over  the  ultraviolet  (UV)  to  far-IR 
spectral  region.  For  devices  such  as  optical  rotators  and  modulators, 
magneto-optical  properties  of  discrete-energy-level  transitions  offer 
significant  advantages  over  bulk  properties  which  are  used  in  conven- 
tional designs. 

In  addition  to  novel  device  applications,  it  has  been  shown  how 
magneto-optical  spectroscopy  can  also  provide  information  about  basic 
parameters  of  a system,  not  easily  obtainable  by  other  means. 
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